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Abstract 


Bacterial disease is a significant issue for larviculture of the Yesso scallop (Patinopecten yessoensis). One source of 
bacteria is the water of larval-rearing tanks. In this study, water bacterial communities were investigated by MiSeq 
sequencing during larval development at fertilized egg W1, trochophora W2, D-shaped larvae W3, umbo larvae W4, 
juvenile scallop W5 stages. Genomic DNA was extracted from the water samples, and a gene segment from the V3-V4 
portion of the 16S rRNA gene was amplified and sequenced using an Illumina Miseq sequencer. A total of 69245 
optimized reads obtained from five samples were identified as 9 phyla and 86 genera. Bacterial communities were 
dominated by the Proteobacteria and Bacteroidetes at phylum level across five larval development stages. At genus level, 
sequences assigned to Pseudoalteromonas, NS3a marine group and Sulfitobacter were dominant at stages of W1, W4 and 
W5. Moreover, reads from Vibrio became dominant at stage W2, followed by Pseudoalteromonas and Sulfitobacter; 
Glaciecola was dominant at stage W3, followed by NS3a marine group and Sulfitobacter. These findings indicated that 


water bacterial communities changed greatly during the trochophora and D-shaped larvae stages of larval development. 


Keywords: Patinopecten yessoensis; Bacterial community; 16S rRNA; High-throughput sequencing 


Introduction 


The Yesso scallop Patinopecten yessoensis has become 
one of the major cultured species in Liaoning and 
Shandong Provinces, China, making an important 
contribution to the socioeconomic development of these 
coastal regions. The larval stage of scallop is susceptible 
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to mass mortality with opportunistic pathogens [1]. The 
successful production of healthy scallop larvae, to a large 
extent, depends on the quality of water including the 
species of bacteria present. Sandaa et al. studied the 
bacterial communities in the Great scallop Pecten 
maximus aquaculture systems using denaturing gradient 
gel electrophoresis (DGGE) [2]. As indicated by DGGE, 
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there are high similarities in the bacterial community 
compositions among different larvae aquaculture systems 
(stagnant, flow-through and a flow- through system with 
seawater treated with ozone) [2]. The dominant 
culturable genera in the water of Chilean scallop 
Argopecten purpuratus larval rearing tanks and the 
survival of larvae are different in the presence and 
absence of antibiotic [3]. 


Massively parallel next generation sequencing 
platforms have now been successfully utilized to 
characterize seawater bacterial communities [4-6]. Sun et 
al. (2016) compared the bacterial communities associated 
with five Yesso scallop larval development stages by 
Mlumina high-throughput sequencing, and found that the 
larval mortality occurred during the transition from the 
D-shaped to umbo larvae and across the entire period of 
umbo larvae stage coincided with the great change of 
bacterial community [7]. Vibrio splendidus and Vibrio 
tasmaniensis are proved to be the pathogens causing 
mortality in Yesso scallop larvae [1]. However, much less 
is known about the bacterial composition in the water of 
Yesso scallop larval-rearing tanks. It is also unclear that 
whether the bacteria in the water are related to the larval 
mortality. The aim of this study was to characterize and 
compare water bacterial communities at different larval 
development stages of the Yesso scallop using Illumina 
high-throughput sequencing, providing the basis for the 
management of the scallop seedling system. 


Materials and Methods 


Sample Collection and Processing 


Water samples were collected from the Yesso scallop 
hatchery in Zhang Zidao Island Group Co. Ltd., Dlian, 
China between March and April 2014. After filtration and 
precipitation, seawater was used in the hatchery. Larvae 
were cultured in plastic tanks containing 100 L seawater 
under static conditions at a density of 5-6 larvae/mL and 
fed Dicrateria inornata 100 mL, 7 - 9x105cells/mL three 
times a day. Low-pressure electrical blowers provided 
aeration via air stones. Water temperature was kept at 15 
+ 1°C and salinity at 32 + 0.5 PSU (Practial Salinity Units). 
Half (50%) of the seawater each tank was replaced every 
day. During larval development stages, three water 
samples of approximately 500 mL each tank were pooled 
and filtered onto 0.22-um filters which were stored frozen 
at -80°C until required for DNA extraction. 


DNA Extraction and PCR Amplification 


Genomic DNA was extracted using a soil DNA 
extraction kit (Shanghai Sangon Biological Engineering 
Technology and Services Company Limited, Shanghai, 
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China) A gene segment (V3-V4 portion) of 16S rRNA gene 
was amplified through PCR by pre-denaturing at 95°C for 
2 min followed by 30 cycles of denaturing at 95°C for 30 s, 
annealing at 55°C for 30 s and extending at 72°C for 45 s, 
and an additional extension cycle at 72°C for 10 min. The 
gene was amplified with primers 338F (5’-barcode- 
ACTCCTACGGGAGGCAGCA-3’) and 806R (5’- 
GGACTACHVGGGTWTCTAAT-3’), in which barcode is an 
eight-base sequence unique to each sample [8-9]. PCR was 
performed in triplicate and the reaction volume was 20 uL 
containing 4 uL of 5 x Fast Pfu buffer, 2 uwL of 2.5 mM 
dNTP (each nucleotide), 0.4 uL of 5 uM primers (each 
direction), 0.4 wL of Fast Pfu polymerase (2.5U/pL), and 
10 ng of template DNA. 


Illumina MiSeg Sequencing 


Amplicons were extracted from 2% agarose gels and 
purified using the Axyprep DNA Gel Extraction Kit 
(Axygen Biosciences, Union City, CA, U.S.) according to the 
manufacturer’s instructions and quantified using 
QuantiFluor™ -ST (Promega, U.S.). Purified amplicons 
were pooled in equimolar and paired-end sequenced (2 x 
300) on an Illumina MiSeq platform following the 
standard protocols (Majorbio Bio-Pharm Technology Co., 
Ltd., Shanghai, China). 


Processing of Sequencing Data 


Raw fastq files were demultiplex and quality-filtered 
using QUME (version 1.17) [10] with the following 
criteria: (i) 300 bp reads truncated at any site over a 50 
bp sliding window were scored in average <20, truncated 
reads shorter than 50 bp; (ii) any mismatched in barcode 
region, 2 nucleotide mismatched in primer and reads 
containing ambiguous characters (iii) only sequences that 
overlap longer than 10 bp assembled, and reads failed to 
be assembled discarded. Operational Taxonomy Units 
(OTUs) were clustered with 97% similarity cutoff using 
UPARSE (version 7.1 http://drive5.com/uparse/) [11] 
and chimeric sequences were identified and removed 
using UCHIME [12]. The most abundant sequence each 
OTU was chosen as the representative. The taxonomy of 
each 16S rRNA gene sequence was analyzed by RDP 
Classifier (http://rdp.cme.msu.edu/) [13] against the silva 
(SSU115)16S rRNA database [14] using confidence 
threshold of 70% [15]. Rarefaction analysis, Coverage (the 
Good’s coverage), Ace (the ACE estimator), Chao (the 
Chao1 estimator), Shannon and Simpson indices were 
created using mothur software (version  v.1.30.1 
http://www.mothur.org/wiki/Schloss_SOP#Alpha_divers 
ity). Venn, mcrobial community barplot, heatmap and 
principal component analysis (PCA) were performed in R 
software. 
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Results 


Illumina Miseq Sequencing Data 


A total of 69245 optimized reads were obtained from 
water samples .Rarefaction curves generated at the OTUs 
level were approached saturation, indicating that the 
sequencing depth was enough to reflect the diversity. 
These reads were clustered into 209 OTUs and the 
number of OTUs of each water sample ranged from 166 to 
188 (Table 1). Good’s coverage estimations revealed that 
99.76% to 99.84% of the species were retrieved (Table 1). 


Our analysis revealed 9 phyla, 17 classes, 33 orders, 50 
families and 86 genera in total. 


Species Richness and Diversity Estimates 


Bacterial taxonomic richness and diversity were 
different among five water samples. As indicated by Chao 
and Ace indices (Table 1), W3 stage harbored the highest 
community richness, followed by W2, W4 and W1 stages, 
while the lowest community richness was found at W5 
stage. Shannon and Simpson indices (Table 1) showed 
that the bacterial diversity increased in order of stages 
W1, W4, W3, W5 and W2. 


Sample OTUs Coverage (%) Ace Chao Shannon Simpson 
W1 170 99.76 189 194 2.67 0.19 
W2 179 99.8 199 196 3.26 0.08 
W3 188 99.84 210 208 2.94 0.12 
W4 179 99.82 194 194 2.8 0.17 
W5 166 99.81 176 173 3.13 0.11 


Table 1: Alpha diversity obtained from five water samples. 


Community Composition 


The shared OTUs were shown via Venn diagram 
(Figure 1). Each stage contained only a minority of special 
OTUs (0~3). Most of OTUs (147~167) were shared 
between stages 


Figure 1: Venn diagram showing OTUs unique to and 
shared by Yesso scallop larval development stages. 
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Nine identified phyla (Proteobacteria, Bacteroidetes, 
Cyanobacteria, Actinobacteria, Planctomycetes, 
Chlamydiae, Verrucomicrobia, Firmicutes and 
Fusobacteria) and two candidate divisions (BD1-5 and 
OD1) were detected among 5 stages, with the bacterial 
16S rDNA sequences dominated by Proteobacteria 
(74.87% ~ 84.53%) and Bacteroidetes (14.49% ~ 
24.43%). The other 7 phyla and candidate division 
OD1were represented by less than 1% of the total 
sequence reads. Sequences from candidate division BD1-5 
and phyla Cyanobacteria, Actinobacteria, Planctomycetes, 
Chlamydiae and Verrucomicrobia were found at all stages. 
Sequences representing candidate division OD1 existed at 
W1, W2 and W4 stages. Sequences from Fusobacteria 
were only detected at W4 stage and those from Firmicutes 
were found at W5 stage alone. 


Analyzing at genus level revealed 86 different 
identified genera in total. Sequence proportions of 12 
genera, NS3a marine group, Roseobacter clade DC-80-3 
lineage, Roseobacter clade AS-21 lineage, SAR11 Surface 1 
and BD1-5 were greater than 1% (Figure 2). Apparent 
shift in Pseudoalteromonas, Glaciecola and Vibrio was 
observed across different larval development stages 
(Figure 2). 


Sequences assigned to Pseudo alteromonas (40.07%, 
37.88%, 25.96%), NS3a marine group (16.62%, 11.27%, 
17.47%) and Sulfitobacter (8.99%, 8.89%, 8.21%) were 
dominant at stages of W1, W4 and W5; followed by those 
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to Alteromonas (3.40%) and Vibrio (2.91%) at stage W1; 
Alteromonas (6.73%) and Phaeobacter (4.72%) at stage 
W4 and Roseobacter clade DC-80-3 lineage (6.20%) and 
Alteromonas (4.55%) at stage W5. However, reads from 
Vibrio (19.33%) became dominant at stage W2, followed 
by those to Pseudoalteromonas(15.31%), Sulfitobacter 
(10.26%), NS3a marine group (8.78%), Aliivibrio (6.10%) 


and Roseobacter clade DC-80-3 lineage (5.40%); while 
sequences belonged to Glaciecola (27.90%) was dominant 
at stage W3, followed by those to NS3a marine group 
(15.40%), Sulfitobacter (7.87%), Polaribacter (5.43%), 
Alteromonas (5.38%) and Roseobacter clade DC-80-3 
lineage (5.32%). 
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Figure 2: Genus abundance (%) of water at Yesso scallop larval development stages. 
Note: sequences abundance less than 1% at a stage was combined into “others”. | 


PCA of Phylogenetic Relatedness 


PCA score plot showed that the first principal 
component (PC1), which explained 74.92 % of the 
variance in the bacterial communities, separated the W3 
stage and other stages (Figure 3). The W4 stage was 
closely related to W1 stage, whereas the W2 stage was 
separate from the other stages along PC2, which 
represented 20.36% of the total variations (Figure 3). 


Overall, the PC1 and PC2 axes explained 95.28% of the 
variations between the different bacterial communities. 
Hierarchically clustered heatmap analysis based on the 
bacterial community profiles at genus level disclosed that 
W1 and W4 stages grouped together firstly, and they then 
clustered with W5, W2 and W3 stages in order (Figure 4). 
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Figure 3: Principal component analysis of bacterial 
communities among five stages at OTU level. 
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Figure 4: Heatmap based on the genus abundance similarity. 
Note: Double hierarchical dendrogram shows the bacterial distribution among 5 stages. 
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The heatmap plot depicts the relative percentage each 
genus (rows) at each stage (columns). The relative 
abundance of bacterial genera is depicted by color 
intensity with the legend appeared at the bottom of the 
figure. Clusters based on the distance of 5 stages along X- 
axis and the bacterial genera along Y-axis are indicated in 
the upper and left of the figure, respectively. 


Discussion 


During different development stages of Yesso scallop 
larvae, the bacterial communities in the water of larval- 
rearing tanks were diverse, especially, at stages W3 and 
W2 (Figure 2-4). One of the reasons may be the change of 
the bacterial compositions from inlet water [2,16]. 
Similarly, the communities in treated clean seawater used 
in the oyster (Crassostrea gigas) hatchery are variable on 
scales of days [16]. Nevertheless, only minor variations in 
the bacterial community composition in DGGE profile 
have been found in the water collected at different days 
during the growth of the Great scallop larvae [2]. 


The bacterial communities observed in water samples 
collected in this study are generally agreement with 
seawater bacterial communities in shellfish hatcheries, 
being dominated by Proteobacteria- (predominantly 
Gammaproteobacteria and Alphaproteobacteria) and 
Bacteroidetes -(predominantly Flavobacteria) affiliated 
sequences. Gammaproteobacteria and/or 
Alphaproteobacteria have been found as dominant 
members of the bacterial populations in the water from 
different Great scallop larvae aquaculture systems [2]. 
Analyses conducted using Gammaproteobacteria specific 
probes have also found significant percentages (9%~ 
48%) of this bacterial group in total counts of water 
samples from Chilean scallop larvae culture tanks [17]. A 
high fraction of the culturable bacteria (60%) are 
assigned to Gammaproteobacteria, as observed by Godoy 
et al., in the water of Chilean scallop larvae reared [3]. The 
bacterial communities in the intake and treated clean 
seawater from a commercial C. gigas hatchery are 
dominated by the Alphaproteobacteria, 
Gammaproteobacteria and Flavobacteria [16]. In addition, 
Gammaproteobacteria, Alphaproteobacteria and CFB 
group (Bacteroidetes) have also been detected from 
seawater in a marine hatchery surrounding waters [18]. 
Of 12 genera (Figure 2) described in this study, 
Pseudoalteromonas and Vibrio have been detected in the 
water from Great scallop larvae aquaculture systems [2], 
Alteromonas is one of the dominant genera in the water of 
Chilean scallop larval-rearing tanks [3], and Glaciecola 
and Polaribacter are common in the seawater used for the 
oyster C. gigas hatchery [16]. Vibrio, Alteromonas, 
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Pseudoalteromonas, Shewanella and Marinomonas can be 
isolated from growing-seawater of Ostrea edulis [19]. 
Sulfitobacter, Phaeobacter, Allivibrio and Cobetia, and 
some OTU unclassified to genus level, such as NS3a 
marine group, Roseobacter clade DC5-80-3 lineage and 
AS-21 lineage have been found from coastal seawater [20- 
24]. Neptuniibacter was first isolated from surface water 
of the eastern Mediterranean Sea [25]. 


In this study, relative abundances of Aliivibrio and 
Vibrio during stages W2 and W3 were higher than those at 
other 3 stages (Figure 2, Figure 4), which might be related 
larval mortality occurred during the transition from D- 
shaped to umbo larvae and at the entire period of stage 
W4. Vibrio splendidus V1 and V. tasmaniensis V2 strains 
have been found to be associated with mass mortality of 
Yesso scallop larvae [1]. The sequences of OTU 102 (V. 
splendidus) at stages W2 and W3 was more abundant than 
that at other 3 stages, which might be the source of larval 
pathogen due to the fact that higher abundance of OTU 
102 from larvae was observed at stages W3 and W4 [7]. V. 
splendidus has also been identified as the pathogen for the 
mass mortalities of Great scallop and Chilean scallop 
larvae [26-27]. A Vibrio anguillarum-related strain 
pathogenic to Chilean scallop larvae has been identified 
[28]. Vibrio pectenicida has been isolated and described 
by Lambert et al. as being the major cause of mortality 
during the production of Great scallop larvae in France. 
Regarding the bacteria in the present study, an 
assessment of the virulence of species in genus Aliivibrio 
in the Yesso scallop larvae is warranted [29]. 


Kesarcodi-Watson et al. has screened probiotics for 
Greenshell™ mussel Perna canaliculus larvae from three 
sources including hatchery seawater [30]. 
Pseudoalteromonas haloplanktis (formerly Alteromonas 
haloplanktis) and Vibrio sp. are able to confer protection 
to Chilean scallop larvae against infection of pathogenic V. 
anguillarum [31-32]. For the larval rearing of Great 
scallop, protection against Vibrio coralliilyticus and V. 
splendidus is provided by prior administration of 
Phaeobacter allaeciensis and Alteromonas macleodii 0444. 
Significant protection with Pseudoalteromonas sp. D411 is 
only provided against V. splendidus [33]. However, 
whether the bacteria in genera of Pseudoalteromonas, 
Vibrio, Phaeobacter and Alteromonas found in this study 
have a probiotic effect in Yesso scallop larvae culture is 
unclear and requires investigation. 


Conclusions 
In summary, water bacterial community compositions 


were diverse across five Yesso scallop larval development 
stages, especially in the middle stages. With our 
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